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We derive the Deep Inelastic anti-quark distribution in a baryon at a low value of using 
the variational principle of Quantum HadronDynamics, an alternative formulation of Quantum 
ChromoDynamics. It is determined by a variational approach generalizing the "valence" quark 
approximation of earlier papers. We find that the "primordial" anti-quarks carry less than a percent 
of the baryon momentum. In the limit of chiral symmetry and Nc <x, we show that the anti-quark 
^\ ' content of the proton vanishes at low Q^. 
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, In previous papers |J1I| one of us has outlined a way of calculating the Deep Inelastic Structure 

" , ' , functions of the baryon from Quantum ChromoDynamics (QCD). After some approximations, ((i) 

. Dimensional Reduction to two dimensions, (ii) Ignoring transverse gluon degrees of freedom) the 

00 ■ 

■ theory reduces to two dimensional QCD, which was transformed into a new form called Quantum 
. HadronDynamics (QHD). In this form the basic dynamical variable is a color invariant quatity 
. M{x,y) = A- ■ [xaix),x^°'iy)] (where X)X^ the annihilation-creation operators of quarks) 

■ which can be thought of as the field operator of a meson. The main advantage of this new point 
^ • of view is that the (iii) semi-classical approximation of QHD corresponds to the large Nc limit of 
^ • QCD, and so is capable of describing non-perturbative phenomena such as the struture of hadrons. 
^ ' The baryon is a topological soliton in this theory and its structure functions (within these ap- 
^ ' proximations) can be determined by a variational principle. In previous papers we made yet an- 

' other approximation, (iv) the assumption of a factorized ansatz for the classical meson variable 

$_( ' M{x,y) = —2tl){x)tp*{y), which corresponds to the valence quark approximation in the parton 



model. 

We have already discussed the consequences of relaxing some of these simplifying assumptions. 
For example the effect of transverse momenta (departure from two dimensionality: relaxing (i)) 
can be studied within perturbative QCD; indeed as emphasized by Altarelli and Parisi, this is the 
physical meaning of the usual DGLAP evolution equations Q for the structure functions. The 
effect of reinstating transverse gluons (relaxing (ii)) is to produce a slightly more involved two 
dimensional field theory, which we will study in a separate paper. (This is important to derive 
the gluon distribution functions of the baryon.) The effect of finite Nc is (in the leading order) to 
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restrict the range of values of the parton momentum [^-^ . In this paper we will study the departure 
from the factorized ansatz for M{x, y); in other words, we will study the departure from the valence 
parton model. This is the same as studying the anti-quark content of the baryon. 

We will see that the probability of finding an anti-quark inside a proton is quite small (< 1%) 
justifying the valence parton approximations made in previous papers. Using a variational ansatz 
we will obtain the anti-quark distribution functions. These can be used as initial data for evolution 
in using the DGLAP equations, which take into account the perturbative corrections. That the 
initial anti-quark content is quite small, is consistent with the phenomenological model of Gliick and 
Reya Q: we now have a theoretical derivation of this picture. However, we expect the initial gluon 
distribution to be non-zero. There are other approaches to studying parton distribution functions, 
see for example Ref. 

Let us begin by summarizing the large Nc limit (which is the classical limit) of QHD. The 
dynamical variable is a complex valued function AI{x, y) of two space-time points (a;, y) lying along 
a null-line. This variable satisfies M*{x,y) — M{y,x) so that we can regard it as the integral kernel 
of a hermitean operator on L^{R). (For technical reasons we assume that this operator is Hilbert- 
Schmidt; i.e., J \M{x,y)\^dxdy < oo.) Moreover, it satisfies the non-linear constraint [e + M]'^ = 1 
where the operator e is the celebrated Hilbert transform operator with — 1 and the integral kernel 
e{x,y) — J sgn (p)e'''''^~"'|^ = j^V-;;:^. This constraint can be understood as a consequence of 
the Pauli principle for fermions as explained in [^. The static solutions of the theory are then the 
minima of the energy functional 

E(M) I f ^-r, ^l r M^" 

4r = -2y^^(^'^)2[^+7j2. ■ 8 

subject to the above constraints. The first term is just the kinetic energy in null co-ordinates 

the second is the potential energy induced by the longitudinal gluon fields. (Recall that the linear 

potential i|a; — j/| is the Fourier transform of the gluon propagator ^ in two space time dimensions). 

The parameter g ~ Aqcd determines the strength of the strong interactions; also, fi^ — — -^^ is 

related to the current quark mass m through a finite renormalization. We will be mainly interested 

in the case m « g which corresponds to the limit of chiral symmetry. It has been shown elsewhere 

^ that subject to the above constraints, the energy E{M) is positive: the constraints being crucial 

for this conclusion. A Lorentz invariant form of the above variational principle is to minimize the 

invariant mass-squared rather than energy. Since the null momentum of a configuration is 
p 



^ + I dxdy\M{x,y)\';-\x-y\ (1) 



■UpMiP,P)t we get 
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2 / P^-^(P'P)^ 



^ ' M{P,P)^^ + 4- I dxdy\Mix,y)\'^\x-y\ 



2 J ''^'^'2p27r 8 J ' 2' 

The quantity B = — | tr M — — ^ J M[x,x)dx is an integer, a topological invariant of the 
configuration as shown in If we reexpress it in terms of the quark fields we can see that this 
is just the baryon number. Thus, a baryon is a topological soliton in this picture. We can get the 
structure functions of the baryon by minimizing the energy functional E{M) subject to the above 
constraints. We have developed a method |^ to solve this problem: a variant of the steepest descent 
method that takes into account the non-linear constraint. However, this method is computationally 
intensive. A method based on a variational ansatz that builds on our previous results on the 



(2) 
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separable ansatz gives as good results while being much simpler computationally. We report on the 
results of this variational approach in this paper. 

The separable ansatz, M{p,q) = —2ip{p)ip*{q) satisfies the constraint if the vector -i/) G L^{R) is 
of unit length and is non-zero only for positive momentum. This is easily verified by substitution. 
In the limit of chiral symmetry m — 0, and Nc ~* ao the exact minimum of the functional is of the separable form, 
with ip{p) ~ e 9. To see this we first derive an integral equation for the extremum of by varying with respect to 



= 0, and can be converted to a 



A/, respecting the constraint. In operator language this equation is 
nonhnear integral equation. By direct computation, the separable ansatz with 'il'{p) ^ e^' can then be verified to be 
an exact solution of this equation, in the limit m = 0. (Details will be given in a longer paper.) Moreover, Ai^ is zero 
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for this solution. Thus the departure from the valence parton picture is determined by the dimensionlcss ratio ^ 
which quantifies chiral symmetry breaking. (As we showed in Ref. Q the leading effect of finite Nc is not to depart 
from the valence parton picture but rather to constrain the range of momenta of the partons.) 

Thus we should expect the 'primordial' anti-quark distribution in the proton to be small: the up and down quarks 
have current quark masses small in comparison to Aqcd; which means that ^ << 1 as well. 

The mathematical advantage of the separable ansatz is that it 'solves' the nonlinear constraint on M: more 
precisely, it replaces it with the condition that is of norm one. In the same spirit, consider the confguration 
M = X^I fc=i CbV'a ^ ip^''- Here we choose tpa to be a set of r orthonormal eigenvectors of the operator e; i.e., 
eipa — Salpa, = ±1- This implics that the operator M is of rank r: the special case of rank one is just the separable 
ansatz above. This ansatz will satisfy the constraint on AI if the r x r matrix ^ is hermitean and satisfies the constraint 
^a^b + [^a + ^c]£,a = 0: a 'mini' version of the constraint on M. Moreover, the baryon number is B = — ^tr AI = —^tr^. 
In the special case of rank one, we have simply ^ — —2. By choosing a large enough value of r this ansatz can produce 
as general a configuration in the phase space as needed: such configurations form a dense subset of the phase space. 

The simplest configuration of baryon number one that departs from the separable ansatz is of rank three. We 
will find that in physically interesting situations, even this departure is very small, so we do not need to consider 
configurations of higher rank. 

By a choice of basis among the ^pa^ we can always bring a rank three configuration of baryon number one to the 
form M = -2ip ® + 2C-{C-[V'- ® "0- - V'+ ® ip+] + - C-liV"- ®ipX + ip+® -0-]} where 0-, V", 0+ are three 
vectors in L'^{R) satisfying eip- = —ip-, ^V' = "07 ~ V-'+JIV'-lP = = II'0+IP = 1) < 07''/'+ >= 0. 

The conditions < ■0-7 >=< 4'-,ip+ >~ are then automatic. The parameter < C- 1 measures the deviation 
from the rank one ansatz and hence, the anti-quark content of the baryon. For example, baryon number is given 
by i? = /o°° |l'0(p)P + C- l0'+(p)P ~ l0-(~p)P I if ■ The wavefunctions 'ip,ip+ both describe quarks and their 
orthogonality can be interpreted as a consequence of the Pauli principle. ^ describes "valence" quarks while is 
the wave function of the "sea" quarks. Since tp_ contributes with a negative sign to the baryon number, it describes 
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anti-quarks. From our previous result we expect C_ to vanish as ^ *■ 0. 

We can substitute this ansatz into the energy E(M) or Al^ and derive integral equations for the minimiza- 
tion. However, in keeping with the spirit of the variational ansatz, we can simplify the problem by assuming 
first some simple functional forms for the functions ip,ip±- The form of the exact solution suggests the choice 
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e '^a for p > and ip-{p) = 'ip{—p) for p < 0. (For other ranges of p 
these functions must vanish.) The parameter Ci is determined by the orthogonahty condition while C, C+ are fixed 
by the normalization conditions. The variational parameter b determines the reference frame. The Lorentz invariant 
quantity A^^ is independent of b. Thus the variational principle will determine a and C- and hence the wavefunctions. 

The rest of the calculation is a straightforward evaluation of the energy integrals and then their minimization. (We 
use the symbolic package Mathematica for some of the computations, most of which can be done analytically. Some 
details are provided in ||^). We have done the calculation and shown that for physically interesting values of ^ 
( A^'^'^ )^ ^ .001), the parameter C_ is quite small. We present the results in the figures which show the small 
effects of deviations from the separable ansatz (i.e. the effects of anti-quarks) and from chiral symmetry. Finally, 
the effect of finite Nc is (in the leading order) to restrict the maximum value of parton momenta. We have already 
studied this correction in the case of the separable ansatz and find it to be small In the case of anti-quarks, 

we establish that they carry less than a percent of baryon momentum in the Nc oo limit and therefore, corrections 
due to finite Nc are less relevant. They will be addressed in a longer paper. 

Thus, we have derived the "primordial" anti-quark distribution function of the proton by a series of approximations 
from QCD. We have an explanation of why it is small in comparison to the valence quark distribution at the low 
initial value of Qq ^ OAGeV^ [^j. The anti-quark distribution is in fact zero in the limit of chiral symmetry and when 
Nc oo, while deviations arc small. This justifies the valence parton approximations made in earlier papers 
It is possible to compare our prediction with experimental data: there is a specific combination of deep inelastic 
structure functions that describes anti-quarks To make a comparison, we need to evolve our distribution from Qg, 
according to the DGLAP equations. However, it is necessary to know the initial gluon distribution in order to solve 
the evolution equations. We will study the gluon distribution in a later paper and subsequently return to this issue. 
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Figure 1: Variational estimates for (a) the exponent 'a' and (b) as a percent. They are plotted as functions of 
nu = 1000 * The exponent a and the anti-quark content (_ go to zero for small current quark masses. 
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Figure 2: (a) Variational estimate for the fraction of fermion momentum of the baryon carried by anti-quarks. It 
is plotted as a percent as a function of nu. The "primordial" anti-quarks carry less than a percent of the portion 



4 



of baryon momentum shared between quarks and anti-quarks, (b) Variational upper-bound on the invariant Mass^, 



( -ajv^i ) ■ of the baryon in the two-dimensional approximation, plotted as a function of nu. In the limit of chiral 



symmetry, we recover the exact exponential solution with 



= 0, C- = and a = 0. 
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Figure 3: (a) The "valence" quark {xb * IV'P)? (b) anti-quark {xb * IV'- P) ^i^d (c) "sea" quark [xb * IV'+P) 
distributions plotted as a function of momentum fraction xb = p at low (~ GAGeV^, see Ref. [^). The 
exponential tails beyond xb = 1 are an artifact of the large- A^c limit. They arc plotted for a small value of current 



quark mass (Sj 



.001) in the reference frame in which the mean baryon momentum, P is 1. The 



fermions are assumed to carry / = 5 the mean baryon momentum. The rest is carried by gluons The node in 

the "sea" quark wavefunction is because it is required to be orthogonal to the "valence" quark wavefunction by the 
Pauli principle. The valence quark distribution shown above (xg* Valence ~ xF^), though calculated in the limit 
Nc —* 00 agrees well with the distribution obtained after taking into account the leading corrections and also with 
experimental measurements of the neutrino structure function xbF^{xb, Q^) when evolved to higher values of 
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